COBE Preprint # 96-04, Submitted to ApJ Letters 



Band Power Spectra in the COBE DMR 4- Year Anisotropy Maps 

G. Hinshaw2'3, A.J. Banday^'^, C.L. Bennett^ K.M. Gorski^'^, A. Kogut^, 
G.F. Smoot'^'S & E.L. Wright^ 

ABSTRACT 

We employ a pixel-based likelihood technique to estimate the angular power 
spectrum of the COBE Differential Microwave Radiometer (DMR) 4- year sky maps. 
The spectrum is consistent with a scale-invariant power-law form with a normalization, 
expressed in terms of the expected quadrupole anisotropy, of Qj.ms-PS\n=i = 18 ± 1.4 
liK, and a best-fit spectral index of 1.2 it 0.3. The normalization is somewhat smaller 
than we concluded from the 2-year data, mainly due to additional Galactic modeling. 
We extend the analysis to investigate the extent to which the "small" quadrupole 
observed in our sky is statistically consistent with a power-law spectrum. The 
most likely quadrupole amplitude is somewhat dependent on the details of Galactic 
foreground subtraction and data selection, ranging between 7 and 10 /iK, but in no 
case is there compelling evidence that the quadrupole is too small to be consistent 
with a power-law spectrum. We conclude with a likelihood analysis of the band power 
amplitude in each of four spectral bands between i = 2 and 40, and find no evidence 
for deviations from a simple power-law spectrum. 
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1. Introduction 

The detection of large angular scale anisotropies in the Cosmic Microwave Background (CMB) 
radiation was first reported by the COBE-DMR experiment in 1992 (Smoot et al. 1992; Bennett 
et al. 1992; Wright et al. 1992; Kogut et al. 1992). The initial detection was based only on the 
first year of flight data. Since that time the DMR Team processed and analyzed the first two years 
of data and found results to be consistent with the first year results (Bennett et al. 1994, Gorski 
et al. 1994, Wright et al. 1994a). We have now processed and analyzed the full 4-years of DMR 
observations: this paper is one of a series describing the results of our analysis. The maps and an 
overview of the scientific results are given in Bennett et al. (1996). 

In this paper we analyze the angular power spectrum of the 4-year DMR maps using a pixel 
based likelihood technique which was first applied to the 2-year data by Tegmark &; Bunn (1995). 

We extend previous work by considering several parameterization of the angular power spectrum. 
The simplest model for large angular scale anisotropy is the power-law model parameterized by 
a normalization, Qrms-PS-, aiid spectral index, n. It is of interest to separate the quadrupole 
anisotropy from the rest of the power spectrum since it is most plausibly contaminated by Galactic 
emission, and, in some models, is predicted to deviate from the simple power-law form. We extend 
our likelihood analysis to a three parameter model in which the quadrupole is fit independent of 
the higher-order power (which is assumed to follow a power-law, see §4 for details). Lastly we 
consider band power estimates in which the anisotropy is assumed to be scale-invariant in each 
of four modestly narrow £ bands, chosen to have roughly comparable sensitivity. The results are 
compared to the power-law fits and indicate that the anisotropy has no significant deviation from 
a power-law form. 



2. Data Selection 

The DMR experiment has produced two independent microwave maps (A and B) at each of 3 
frequencies (31.5, 53 and 90 GHz). The results presented here are based on linear combinations of 
all 6 channel maps. The combination coefficients are dictated by the sensitivities of the individual 
channels and on considerations of Galactic foreground removal. Kogut ct al. (1996a, 1996b) have 
revisited the question of Galactic emission in the 2- and 4- year DMR data at high latitudes, and 
conclude that there is statistically significant evidence for a weak Galactic signal at all three 
frequencies, even at latitudes |6| > 20°. Thus, for the 4-year analysis, we take a more aggressive 
approach to Galactic foreground removal than we have previously: first, we extend the Galactic 
plane cut of 20° with additional cuts, guided by the COBE-DlY(SE 140 jJLva map (Bennett et 
al. 1996). The number of pixels surviving the cut is 3881 in Galactic coordinates. Second, we 
model and remove residual high-latitude Galactic emission in two complementary ways, described 
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below. In all, we analyze three separate maps in this paper: the first map is a weighted average 
of all six DMR channel maps with no residual Galaxy emission subtracted. We denote this map 
"31+53+90". The second is the same weighted average map as the first with best-fit Galaxy 
template maps subtracted from each channel prior to averaging (Kogut et al. 1996b). We denote 
this the "Correlation" model map. The third is a linear combination of all six channels with 
coefficients designed to maximize sensitivity subject to the constraint that any signal with a 
free-free frequency spectrum {Pff = —2.15) cancels. This map also has best-fit synchrotron and 
dust emission templates subtracted prior to averaging (Kogut et al. 1996b). Wc denote this the 
"Combination" model map. The specific coefficients used to construct these maps are given in 
Table 1 of Hinshaw et al. (1996). In all the analyses below we use the maps pixelized in Galactic 
coordinates. Gorski et al. (1996) and Banday et al. (1996) have analyzed both the Galactic and 
ecliptic maps in detail. Since our results for the Galactic maps are in agreement with theirs, where 
comparable, we defer to those papers for a comparison of the Galactic and ecliptic maps. 



3. Method 



Most cosmological models make predictions for the mean angular power spectrum of 
CMB anisotropies, the coefficients, Q. For a rotationally invariant theory, the Q specify the 
expected variance in each spherical harmonic mode in a Fourier expansion of the sky temperature 
T{0, (j)) = Y.e,m aemYem{0, 4>) with {aimO-i'm') = ^mm'- For a given power spectrum, Q, the 

implied covariance between map pixels i and j is given by 

M,,- = {TiTj) = ^ (2£ + 1) W! Q Piim ■ rij) (1) 
e 

where Tj is the temperature in pixel z of a map, the angled brackets denote a universal ensemble 

average, is the experimental window function that includes the effects of beam smoothing and 
finite pixel size, Ci is the power spectrum, Pi{ni ■ rij) is the Legendre polynomial of order i, and lii 
is the unit vector towards the center of pixel i. For Gaussian fluctuations, the covariance matrix 
fully specifies the statistics of the temperature fluctuations. The probability of observing a map 
with pixel temperatures T, given a power spectrum Ci, is 

df e-i3"-M(C<)-'-J' 

where N is the number of pixels in the map. Assuming a uniform prior distribution of cosmological 
model parameters, the probability of a power spectrum C^, given a map T, is then 

lfT.jv/-l(C,).f 

£(Q|r) oc ^ -■ (3) 

^ ' ^ VdetM Q ^ ^ 
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In the following section, we evaluate the above likelihood function using three different 
parameterizations of the power spectrum, C^. To test the effects of data selection and Galaxy 
modeling, we analyze three separate DMR maps, as specified in §2. To make the analysis 
computationally efficient, we have degraded the maps by one step in pixel resolution (to index level 
5) for which there are 1536 pixels in the full sky, and 954 pixels surviving the extended Galaxy 
cut. We account for the effects of smoothing due to pixelization by including a term in the window 
function: Wi = G^Fi. The are the Legendre coefficients of the DMR beam pattern, tabulated 
by Wright et al. (1994b). The i*) are the Legendre coefficients for a circular top-hat function with 
area equal to the pixel area. The coefficients for index level 5 pixels are available on request. 

We ignore the contribution of the monopole and dipole moments in the maps since these 
modes are either unconstrained by the data (the monopole), or are dominated by local effects (the 
dipole). In principle, this is achieved by integrating the likelihood over the modes Cq and Ci; in 
practice, we set these terms to large positive constants when evaluating the covariance matrix. 
We have found that setting Cq = Ci = 10® /xK^ renders the likelihood insensitive to monopole 
and dipole moments of several hundred /iK, without compromising the inversion of the covariance 
matrices. (Note that the analyzed maps have approximately zero mean, by construction, and 
approximately zero dipole since an estimate of the CMB dipole is removed during the raw data 
processing.) We assume the noise in the sky maps is uncorrelated from pixel to pixel (Lineweaver 
et al. 1994), which adds a diagonal contribution to the pixel covariance matrix in Equation |l[ 
Tegmark &; Bunn (1995) have shown the assumption of uncorrelated noise to be an excellent 
approximation for this application. The noise per pixel is derived from the noise per observation, 
given in Table 1 of Bennett et al. (1996), and the number of observations per pixel. 



4. Results 



We consider three parameterizations of the angular power spectrum. First, we adopt the 
power-law model, parameterized by the amplitude of the mean quadrupole anisotropy, Qrms-PS-, 
and the power-law spectral index n. Specifically (Bond &: Efstathiou, 1987) 

r an ^\-(A^K^n^ r(£ + (n - i)/2)r((9 - n)/2) 

Ce = CeiQrms-PS, n) = (4vr/5)Q__p5 r(^ + (5 - n)/2)r((3 + n)/2) 

This model is extended to study the quadrupole anisotropy by parameterizing the power at £ = 2 
separately 

\ Ce{Qrnis-PS,n) ^ > 3 

The most-likely value of C2 that results from this model is closely related to the quadrupole 
anisotropy observed in our sky, which we denote Qrms] the precise connection is discussed below. 
Note also that the power-law parameters, Qrms-PS and n, inferred from this model are essentially 
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the same as those derived from marginahzing over C2 since they are only weakly coupled to C2. 
Lastly, we study a model in which the spectrum is taken to be scale-invariant in each of four 
relatively narrow i bands, and let the amplitude in each be a free parameter 



Note that within each band, the amplitude parameters QaS correspond to Qjiat as defined by 
Scott, Silk & White (1995). The spectral band widths were chosen to give roughly equal sensitivity 
in each band except the highest which suffers loss of signal due to the 7° beam width. The model 
is designed to probe for deviations from a power-law while maintaining computational feasibility. 

The fits to power-law spectra, including the quadrupole in the analysis, are summarized in 
Table |l]. The results are generally consistent with the 2-year data. The overall normalization is 
slightly smaller due to the additional Galactic cutting and modeling. The most-likely spectral 
index is slightly greater than unity, while the quadrupole normalization for a scale-invariant 
spectrum ranges from 17.2 to 18.4 fiK, depending on Galactic model. For comparison, the 
scale-invariant normalization derived from the weighted average map using a straight 20° cut with 
no additional Galactic modeling is 20.1 ^K, 1.6 fiK higher than we obtain with the extended cut, 
and comparable to the normalization quoted by Gorski et al. (1994) for the 2-year data (using a 
straight cut). In assessing the results obtained from the three maps, we note that the DIRBE 140 
/im map appears to trace the bulk of the free- free emission seen by DMR (Kogut et al. 1996b), 
and since the Correlation map is more sensitive than the Combination map, we give it more weight 
in our conclusions. Taken together, the results in Table 1 are consistent with a spectral index of 
1.2 ± 0.3 and a scale-invariant quadrupole normalization of 18 it 1.4 fiK. 

The fits to power-law spectra with the quadrupole parameterized independently are 
summarized in Table ^. The first three columns summarize the power-law portion of the spectrum 
in the same format as Table ||. These results are based on slicing the 3-dimensional likelihood at 
the maximum likelihood value for C2, but the results are only weakly dependent on C2 and thus 
are effectively equivalent to standard power-law fits that ignore the quadrupole. The last column 
of Table | gives the 68% confidence interval for C2 expressed in terms of Qrms = 
The mode gives a self-consistent estimate of the quadrupole moment observed in our sky, while 
the confidence range accounts for both instrument noise and cosmic variance. A complimentary 
approach to analyzing the quadrupole, based on fitting and squaring a2m coefficients (Kogut et al. 
1996b) gives consistent results, after accounting for the bias introduced by uncertainties in the 
a2m- This approach demonstrates the importance of modeling the residual galactic foreground 
emission since it contributes significantly to the quadrupole emission. As with previous analyses, 
the observed quadrupole is smaller than that expected from the power-law fits: the most-likely 
amplitude ranges from 6.9 to 10.0 depending on Galaxy model. Figure ^ shows the full 
likelihood for Qrms for each map analyzed. It is important to stress that while the quadrupole in 



' (247r/5)Q2 
(24^5) Q2 
(247r/5)Q2 
(24^/5) Q2 



2<i<5 
6<e<10 
11 < ^ < 20 
21 < ^ < 40 



(6) 
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our sky is most likely ~10 ^K, the cosmic variance combined with experimental uncertainties are 
so large that its value is easily consistent with a power-law model of anisotropy. For example, the 
likelihood for Qrms derived from the Correlation map implies there is a 22% chance that Qrms 
exceeds 18 ;uK, the value favored in a scale-invariant power-law model. 

The band power fits are summarized in Table ^, and are plotted in Figure |2|. The vertical 
uncertainties in the figure indicate the extent of the 68% confidence interval in each band power 
parameter when the other 3 are fixed at their maximum likelihood value. These uncertainties 
include both instrument noise and cosmic variance. The horizontal errors bars represent the 
extent of each band, as defined in Equation ^. The covariance between bands, which arises from 
the Galaxy cut and from non-uniform sky coverage, is quite small: roughly 10% of the parameter 
variance for neighboring bands, and less for non-neighboring bands. Note that each of the three 
lowest bands have consistently significant detections of power, while, in all cases, the highest 
band, from 21 to 40, does not. Thus, we only plot 95% confidence upper limits for this band. To 
compare the band power fits to the power-law fits we have also plotted the 68% confidence locus 
of acceptable power-law models in Figure ^. More precisely, the dashed white line in the figure is 
the mean power spectrum for the most-likely power-law model, while the grey band represents 
the locus of mean power spectra within the 68% confidence region in the {Qrms-PSif^) plane. 
The general agreement between the power-law model and the band power model is an indication 
that there are no significant wide band deviations (with ~ a few) from a simple power-law 
in the low-£ anisotropy spectrum. The band power amplitude in the highest d. band we probe is 
consistently low, but this estimate is rather sensitive to the details of the beam and pixelization 
filters, and to the level of noise in the maps, so the uncertainty attached to this estimate is quite 
large. To date, two other experiments have measured anisotropy on angular scales probed by the 
DMR: FIRS (Ganga et al. 1994) and Tenerife (Hancock et al. 1994). Both experiments report 
significant detections of anisotropy: the FIRS team quotes Qfiat = 19 ± 5 //K for ^^30, while the 
Tenerife team quotes Qfiat = 26 it 6 ;uK for 13$ £^30, both of which are consistent with DMR. 



5. Conclusions 

We have estimated various parameterizations of the angular power spectrum in the COBE- 
DMR 4- year sky maps. We find the results to be generally consistent with the first and second year 
results. The data are consistent with a scale-invariant spectrum with a quadrupole normalization 
of 18 lb 1.4 jiK, and a best-fit spectral index of 1.2 it 0.3. The quadrupole anisotropy is somewhat 
smaller than the best-fit power-law spectrum would prefer, but the discrepancy is not statistically 
significant when we take account of Galactic modeling uncertainties, instrument noise, and cosmic 
variance. We have further analyzed the spectrum in each of four I bands and find no evidence for 
significant, wide-band deviations from a simple power-law form. 

We gratefully acknowledge the many people who made this paper possible: the NASA Office 
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analyze the data. We also thank Charley Lineweaver for useful comments on our manuscript. 
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Table 1. Power-law Spectral Parameters 



Map Qrms-PS° Qrms-PS\n=l 

iliK) (/xK) 



31+53+90'^ l-25t[i:i^ 15.4liJ 18.4^1:^ 

Correlate ^■'^'^-o^S '^^■"^-2! 17-8li| 

Combination^ l-00l[ii^ 17.2l|:^ 17.2+^;^ 



'^Mode and ±68% confidence interval for the projection of the 2- 
dimensional likelihood C{Qrms-PSin) on to n. 

^Mode and ±68% confidence interval for the projection of the 2- 

dimensional likelihood C{Qrms-PSin') on to Qrms-PS- 

^Mode and ±68% confidence interval for the slice of the 2- 

dimensional likelihood C{Qrms-PS,n) at n = 1. 

'^Linear combination coefficients for the maps analyzed here are 
given in Table 1 of Hinshaw et al. (1996). 
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Table 2. Quadrupole + Power-law Spectral Parameters 



Map 




Qrms—PS 


Q rrns — P S\n=l 


O 

^ rrns 








(mK) 


(mK) 


31+53+90^ 


1 n9+°-29 

-'-•'-'^-0.30 




18.7tli 


6.9l|-^ 


Correlation^ 


1 Aq+0.29 
-'-•'-'^-0.31 


17.011:^ 


18.lt\i 


lo.ollt 


Combination^ 


n (-7+0.44 
^•'^'-0.49 


23.011^ 


17.9111 


7.610 



^Mode and ±68% confidence interval for the projection of the 2- 
dimensional likelihood C{Qrms-PS-i'n) on to n. C{Qrms-PS^'n) is 
the 3-dimensional likelihood C{Qrms-iQrms-PSi''^) evaluated at the 
maximum likelihood value of Qrms- 

'^Mode and ±68% confidence interval for the projection of the 2- 

dimensional likelihood C{Qrms-PS-, n) on to Qrms-PS- 

^Mode and ±68% confidence interval for the slice of the 2- 

dimensional likelihood C{Qrms~PS^''^) at n = 1. 

"^Mode and ±68% confidence interval for the slice of the 
3-dimensional likelihood jC(Qrms, Qrms-PS, n) at the maximum 
likelihood values of Qrms-PS and n. Similar results are obtained 
by marginalizing over Q 

rms— PS and n. 

'^Linear combination coefficients for the maps analyzed here are 
given in Table 1 of Hinshaw et al. (1996). 
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Table 3. Band Power Spectral Parameters^ 



Map Band 

2<i<5 6<£<10 11<^< 20 21 <^<40 



31+53+90*^ IS.etii 16.711:^ 20.311:? 1-01l6^ 

Correlation^ IS.OtH 15.9li| IS-Sllo O-S^as^ 
Combination*^ ILbtfl^ 17.2li| 17.2lt^ O.llg^i^ 



^Mode and ±68% confidence interval for the band power 
amplitudes, expressed in terms of Qfiat- Qfiat is the quadrupole 
normalization expected for a scale-invariant power-law spectrum 
within the specified range of i. The units are ^uK. 

^Linear combination coefficients for the maps analyzed here are 
given in Table 1 of Hinshaw et al. (1996). 
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10 20 30 

^rms = (mK) 

Fig. 1. — The likelihood function for the mean quadrupole moment observed in our sky for the 
three maps defined in §2. The curves include cosmic variance and instrument noise. In all cases 
the most likely quadrupole is smaller than that favored by the power-law fits to the full data, but 
the likelihoods are all sufficiently broad to encompass the case Qrms = Qrms-PS- The effect of 
Galactic modeling on Qrms is relatively modest, but it does have a significant effect on the phase 
of the quadrupole, particularly the coefficient 020 (Kogut et al. 1996) 
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Fig. 2. — Power spectral estimates for three maps defined in §2. The points with error bars give the 
most hkely band power amphtude within each t band, as indicated by the horizontal error bars, 
under the assumption that the power spectrum is scale-invariant within each band. The fourth 
band, 21 < ^ < 40, is plotted as a 95% CL upper limit, since there is no significant detection of 
power in this band. The vertical errors include both noise and cosmic variance. The shaded region 
indicates the locus of the mean of acceptable power-law models, i.e., those models within the 68% 
confidence region in the {Qrms-PSj ^) plane, as determined from the power-law fits to each map. 



